Spin Qubits in Multi-Electron Quantum Dots 
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We study the effect of mesoscopic fluctuations on the magnitude of errors that can occur in 
exchange operations on quantum dot spin-qubits. Mid-size double quantum dots, with an odd 
number of electrons in the range of a few tens in each dot, are investigated through the constant 
interaction model using realistic parameters. It is found that the constraint of having short pulses 
and small errors implies keeping accurate control, at the few percent level, of several electrode 
voltages. In practice, the number of independent parameters per dot that one should tune depends 
on the configuration and ranges from one to four. 
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I. INTRODUCTION 

Since the discovery that quantum algorithms can solve 
certain computational problems much more efficiently 
than classical onesji^ attention has been devoted to the 
physical implementation of quantum computation (QC). 
Among the many proposals, there are those based on 
the spin of electrons in laterally confined quantum dots 
(QD), 3 which may have great potential for scalability and 
integration with current technologies. For any successful 
proposal, one must be able to perform single- and double- 
qubit operations much faster than the decoherence time. 
In fact, all logical operations required for QC can be re- 
alized if these elementary operations are sufficiently error 
free^ 

Single qubit operations involving a single QD will likely 
require precise engineering of the underlying material or 
control over local magnetic fields^ both have yet to be 
achieved in practice. Two-qubit operations, in contrast, 
are already within experimental reach. They can be per- 
formed by sending electrical pulses to modulate the po- 
tential barrier between adjacent QDs. That permits di- 
rect control over the effective, Heisenberg-like, exchange 
interaction between the qubit spins, which is created by 
the overlap between the electronic wave-functions of the 
QDs£ These operations are important elements in form- 
ing a basic two-qubit gate such as the controlled-no^ and 
in the propagation of quantum information through QD 
arrays^ In fact, using three QDs instead of just one to 
form a logical qubit would allow one to perform all logical 
operations entirely based on the exchange interaction^ 
Thus, exchange operations will likely play a major role in 
the realization of QD qubits. A quantitative understand- 
ing of errors that occur during an exchange operation will 
help in designing optimal systems. 

The first proposal for a QD spin qubil£ relied on hav- 
ing a single electron in a very small laterally confined 
QD. One advantage of such a system is that the Hilbcrt 
space is nominally two-dimensional. Leakage from the 
computational space involves energies of order either the 
charging energy or the single-particle excitation energy, 
both of which are quite large in practice (~ 1 meV 



~ 10 K). Working adiabatically - such that the inverse of 
the switching time is much less than the excitation energy 
- assures minimal leakage. The large excitation energy 
implies that pulses of tens of picoseconds would be both 
well within the adiabatic regime and below the dephas- 
ing time (which is typically in the nanosecond range 
since orbital degrees of freedom are involved). However, 
in practice, it is difficult to fabricate very small tunable 
devices^ Moreover, one-electron QDs may offer little pos- 
sibility of gate tuning due to their rather featureless wave 
functions. 

Alternatively, a qubit could be formed by the top 
most "valence" electron in a QD with an odd number 
of electrons^ In this case, electrons filling the lower en- 
ergy states should comprise an inert shell, leaving as the 
only relevant degree of freedom the spin orientation of 
the valence electron. Large QDs with 100-1000 elec- 
trons, while much simpler to fabricate than single elec- 
tron QDs, are unsuitable because the excitation energy 
is small (~ 50 /LteV ~ 0.6 K), leading to leakage or ex- 
cessively slow exchange operations. On the other hand, 
mid-size QDs, with 10-40 electrons, are sufficiently small 
to have substantial excitation energies, yet both reason- 
able to fabricate and tunable through plunger electrodes. 
For these dots, a careful analysis of errors is necessary. 

Perhaps the best example of an exchange operation 
is the swap of the spin states of the two qubits. For 
instance, it causes up-down spins to evolve to down-up. 
Maximum entanglement between qubits occurs when half 
of a swap pulse takes place - a square-root-of-swap op- 
eration. Several authors have treated the problem of 
swap errors in QD systems ii2*ii*i2iiii4 A primary con- 
cern was the occurrence of double occupancy (when both 
valence electrons move into the same QD) during and 
after the swap. However, no study so far has consid- 
ered another intrinsic characteristic of electronic states in 
multi-electron QDs, namely, their marked dependence on 
external perturbations such as electrode voltage or mag- 
netic field. This sensitivity gives rise to strong samplc- 
to-sample fluctuations arising from the phase-coherent 
orbital motion, 15 These features can make the precise 
control of energy levels, wave functions, and inter-dot 
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couplings a difficult task. 

In this work we study errors and error rates that can 
take place during the exchange operation of two spin 
qubits based in multi-electron QDs. We consider realis- 
tic situations by taking into account an extra orbital level 
and fluctuations in level positions and coupling matrix el- 
ements. These lead to deviations from a pre-established 
optimal swap operation point, especially when a single- 
particle level falls too close to the valence electron level. 
Reasons for such fluctuations can be, for instance, (i) 
the lack of a sufficient number of tuning parameters (i.e., 
plunger electrodes) , or (ii) the cross-talk between the tun- 
ing electrodes. Our results set bounds on the amount of 
acceptable detuning for mid-size QD qubits. 

This paper is organized as follows. In Sec. [H] we in- 
troduce and justify the model Hamiltonian. The states 
involved in the exchange operation are presented in Sec. 
IIIII where we also discuss the pulses and the parame- 
ters involved in the exchange operations. In Sec. IIVI we 
present the results of our numerical simulations. We also 
discuss the impact of mesoscopic effects on errors and 
put our analysis in the context of actual experiments. 
Finally, in Sec. 0we draw our conclusions. 



II. MODEL SYSTEM 

We begin by assuming that the double QD system can 
be described by the Hamiltoniarii& 

H = H A + H B +H AB , (1) 

where 

H a = ^ £ j Ua + ~Y X! Ua ^ a ( X! n »M' - l) ) ( 2 ) 
j.a j,a k.cr' 

a = A, B, and 

Hab = J2 a \° hk ° + h - c -) ' ( 3 ) 

Here, TiAja = a ]a a j<? anc ^ n B,ka = b\ a bk<j are the number 
operators for the single-particle states in the QDs (named 
A and B), e" denotes the single-particle energy levels, tjk 
are the tunneling amplitudes between the dots, and U a is 
the charging energy (er =| , | and j, k run over the single- 
particle states). Typically for mid- to large-size QDs, the 
charging energy is larger than the mean level spacing. 

In the literature of Coulomb blockade phenomena in 
closed QDs, the Hamiltonian in Eq. J3J is known as 
the constant interaction model. It provides an ex- 
cellent description of many-electron QDs, being sup- 
ported by both microscopic calculations and experimen- 
tal datai 15 i 16 ' 17 i 18 ' 19 The reasoning behind its success can 
be understood from two observations. First, mid- to 
large-size QDs, with more than ten electrons, behave very 
much like conventional disordered metals in the diffusive 



regime. Wavelengths are sufficiently small to resolve ir- 
regularities in the confining and background potentials, 
leading to classical chaos and the absence of shell effects 
in the energy spectrum. In this case, the single-particle 
states obey the statistics of random matrices, showing 
complex interference patterns and resembling a random 
superposition of plane waves. This is in contrast with 
the case of small, circularly symmetric, few-electron QDs, 
where shell effects are pronounced 

Second, for realistic electron densities, the QD linear 
size is larger than the screening length of the Coulomb 
interactions. In the presence of random plane waves, the 
screened interaction can then be broken up into a leading 
electrostatic contribution characterized by the QD capac- 
itance plus weak inter-particle residual interactionsiiSiiS 
This description becomes more accurate as the number 
of electrons gets larger since the residual interactions be- 
come weaker. The electron bunching is reduced as wave 
functions become more uniformly extended over the QD. 
Also, the increase in the number of oscillations in the 
wave functions leads to a self-averaging of the residual 
interactions. In this limit, one arrives at the so-called 
"universal Hamiltonian" for QDs, containing only single- 
particle levels, the charging energy, and a mean-field ex- 
change termii£*iii This Hamiltonian can be derived ex- 
plicitly via a random-phase approximation treatment of 
the Coulomb interaction and the use of random-matrix 
wave functions>iSiiS 

According to these arguments, interaction effects be- 
yond the charging energy term are omitted in Eq. © . In 
addition, the intra-dot exchange interaction, which tends 
to spin polarize the QD, is also neglected. The reason for 
that is the following. One can show that the intra-dot 
exchange term only affect states where there is double 
occupancy of a level. Thus, the exchange interaction 
constant always appears side-by-side with the charging 
energy. But in multi-electron dots, the exchange energy 
(which is at most of order of the mean level separation) is 
much smaller than the charging energy. Thus, intra-dot 
exchange effects are strongly suppressed by the charg- 
ing energy. We have verified that their inclusion does 
not modify appreciably our final results. We expect the 
exchange interaction to become important for two-qubit 
operations only in the case of small QDs with only a 
few electrons, when all energy scales (including the mean 
level spacing) are of the same order. 

Thus, the simple picture where single-particle states 
are filled according to the Pauli principle up to the top 
most level is an appropriate description of multi-electron 
J l ;s> i5 4 i6 i i7 i i8 i i9 j n orc i er t define the spin-i qubits, 

both QDs should contain an odd number of electrons 
(say, 2Na~ 1 and 2Ab — 1). The QD spin properties are 
then dictated by the lone, valence electron occupying the 
highest level. The remaining electrons form an inert core, 
provided that operations are kept sufficiently slow so as 
not to cause particle-hole excitations to other levels. 

Experimentally, the two-qubit exchange operations 
also require the capability of isolating the QDs from each 
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FIG. 1: Schematic disposition of energy levels of a system 
of two QD spin qubits (only levels close to the top occupied 
state are shown) . The dashed lines indicate the most probable 
transitions that can occur during the exchange operation. 

other, so that a direct product state can be prepared, 
such as 



\N A J) A (g>\N B ,i)i 



(4) 



where the kets represent only the spin of the valence elec- 
tron on each QD. 



III. ERRORS IN EXCHANGE OPERATIONS 

We focus our study on errors that appear after a full 
swap operation (which should result in no entanglement). 
Although it could in principle seem more sensible to look 
at the square-root-of-swap operation (which creates en- 
tanglement and is therefore a building block of logical 
gates), error magnitudes for the latter are straightfor- 
wardly related to those of the full swap operation, as we 
will show. We leave the discussion of the square root of 
swap to Sec. IIVI 

The ideal full swap operation exchanges the valence 
electrons of the QD system. For instance, it takes the 
product state \i) into 



I/) = Usw \i) = \Na, I) a ® \Nb, T)j 



(•5) 



Physically, the full swap can be implemented by start- 
ing with isolated QDs, turning on the inter-dot coupling 
for a time T (the pulse duration), and then turning it 
off, isolating the QDs again. For weakly coupled QDs 
< U,5e), one finds T « (tt/4) U/\t\ 2 , where t and U 
here represent typical values for the coupling matrix ele- 
ment and the charging energy, respectively (throughout 
we assume h = 1). To quantify the amount of error that 
takes place during the operation, we use the probability 
of not reaching | /) asymptotically, namely, 



e = l-K/KK+oo))| s 



(0) 



We solved numerically the time-dependent Schrodinger 
equation that derives from Eq. for the particular but 
nevertheless realistic case shown in Fig. We assumed 
that voltage tuning allows one to place the top most elec- 
tron of QD A into an isolated single-particle state of en- 
ergy eff A — aligned with the energy of the top most 
electron in QD B, ef [B = 0. However, limited tuning abil- 
ity leaves an adjacent empty state close in energy in QD 



B: £% B+ i = A. Therefore, while we can approximately 
neglect all levels but one in QD A, for QD B we needed 
to take two levels into account, having hopping matrix 
elements denoted by t\ — tjq A N B an d t% = tN A ,N B +i~* 
To facilitate the analysis, we assumed that the dots have 
the same capacitance, C, so that Ua = Ub = U = e 2 /C. 

The Hamiltonian of Eq. conserves total spin. As- 
suming that filled inner levels in both QDs are inert 
(forming the "vacuum" state |0)), we can span the S z — 
Hilbert subspace with nine two-electron basis states. Ac- 
cording to their transformation properties, they can be 
divided into "singlet" 



(7) 



\ D i) = 5 L+;-i,i 6 k+;-i,Tl°>> 

and "triplet" 



IT,) 



\D,) = 



AJ i^N B +l-l,l a N A ] + b N B +l-l,] a N A l) I )' 



V2 



~M ( b N B + lA b N B T + b N B + l,1 b N B l) 1°)' ( 8 ) 
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classes, with I = 1, 2. 

The final states that correspond to an error have ei- 
ther double occupancy (|-Dfc), k = 1, . . . , 5), or an elec- 
tron in the (N B + l)-level of QD B (\S 2 ) and |T 2 )). 22 
In addition, a return to the initial state is also consid- 
ered an error. It is worth noticing the difference between 
our treatment of the problem and that of Ref. EI In 
our case, errors come mainly from either ending in the 
excited single-particle state after the operation is over 
(i.e., states IS2) and IT2)), or from "no-go" defective op- 
erations. In Ref. Il3l errors come from having double 
occupancy in the final state. Double occupancy errors 
can be exponentially suppressed by adiabatically switch- 
ing the pulse on and off on time scales larger than the 
inverse charging energyAii*H Making pulses adiabatic 
on the time scale of the inverse mean level spacing for 
multi-electron quantum dots is more challenging, espe- 
cially because the spacings fluctuate strongly both from 
quantum dot to quantum dot and upon variation of any 
external parameter (mesoscopic fluctuations). Therefore, 
multi-electron quantum dots require extra tunability to 
get around such problems. 

Very small errors, below 10 -6 -10 -4 , can, in principal, 
be fixed by the use of error correction algorithms. 23 The 
pulses, therefore, should be sufficiently adiabatic for er- 
rors to remain below this threshold. We adopted the 
following pulse shape: 

, . If t + T/2 , t-T/2\ 
v(t) - - (^tanh —^L- - tanh —^J- J , (9) 
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FIG. 2: Full swap error as a function of upper level detuning in 
quantum dot B. The pulse width is optimized for A = 1, r = 
6, and tip = 0.2. Results for different charging energies are 
shown. Interference between different quantum mechanical 
paths in the device causes a sharp minimum. 

where r is the switching time. The pulse will remain 
both well-defined and adiabatic provided that T ^> t ^> 
max{A _1 , U^ 1 }. Notice that this pulse is equivalent to 
that adopted in Ref. |r| up to exponential accuracy, 
0{e- T / T ), with T > r. There is no particular reason 
to believe that either performs better than the other; our 
choice was dictated by technical convenience. 



IV. RESULTS 

We used a standard numerical method, the so-called 
Richardson extrapolation^ to solve the Schrodinger 
equation for \ip(t)}. The first step in our analysis was 
to find the optimal value of T which minimized the full 
swap error, as defined in Eq. JJJJ, for a given set of pa- 
rameters U, t\, and t (we used A = 1 and took t 2 =t\). 
The second step was to study how this minimal error 
depends on r. There is actually an optimal interval for 
r, since small switching times spoil adiabaticity, while 
large ones compromise the pulse shape (when T is rela- 
tively short). Empirically, we find that errors related to 
switching times become negligible once r reaches values 
of about r = 4max{A~ 1 , t/ -1 }, provided that r CT. 
In what follows, we fix r > tq. 



A. Mesoscopic Effects 

Figure [21 shows the full swap error as a function of A 
when T is fixed to its optimal value for A = 1. Such a 
situation would arise experimentally if the pulse is op- 
timized for a certain configuration, but a fluctuation in 
level spacing occurs. Notice the sharp increase in error 
as A decreases. While increasing r reduces this error (by 
making the switching more adiabatic), very small level 
spacings would be problematic, since r can not be larger 
than T without sacrificing pulse shape and effectiveness. 
In order to make space for an adiabatic switching time 
for small A, one would also have to increase pulse du- 
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FIG. 3: Full swap error as a function of detuning in the cou- 
pling constant ti. The parameters used in the pulse width 
optimization are the same as in Fig. [5] Results for different 
charging energies are shown. 

ration. This is clear in the case of U = 1 (see Fig. 
Even moderate couplings, ti f2 = 0.2, lead to larger er- 
rors, which can then be suppressed by decreasing ti$ by 
a factor of two; however, that causes a fourfold increase 
in pulse width which may be problematic in terms of de- 
coherence. 

The dependence of errors on fluctuations in the cou- 
pling amplitude t% is shown in Fig. [21 Again, the pulse 
duration used is the optimal value obtained when t 2 — 
t\ = 0.2. As expected, the error grows as t 2 increases. 
Errors related to large values of t 2 can also be minimized 
by increasing the switching time, but the same issues 
raised above appear. 

Figure^lpresents the error for two situations involving 
pulses with duration of about T/2, corresponding to the 
square-root-of-swap operation. The cases shown are: (i) 
one, and (ii) two consecutive square-root-of-swap pulses. 
For comparison, the curve corresponding to a full swap 
pulse is also shown. The error for the square-root-of-swap 
operation is given by Eq. 10 with |/) replaced by 



I/') 



f 



\Si) 



1 



\Ti) 



(10) 



One can observe from Fig. 01 that error rates are nearly 
the same after a full swap operation and after two consec- 
utive square root of swap operations. This insensitivity 
of the error to the pulse duration led us to concentrate 
our effort on the full swap operations. 

In order to establish an upper bound for QD tuning 
accuracy, we have performed simulations where both A 
and t 2 were allowed to vary. The spacing between the 
levels in QD B was taken from a Gaussian distribution 
centered at A = 1, with standard deviation a a ( € n b = ^ 
was kept fixed). For the coupling amplitude, we gen- 
erated Gaussian distributed level widths T 2 = 2irt\/K, 
with average T 2 — 2itt\/K and standard deviation or 2 . 
The pulses had their widths optimized for the typical case 
where A = A = 1JJ = 4, t 1)2 = 0.2, and r = 6. For fixed 
values of A, a a, r 2 , and or 2 , we generated 10,000 real- 
izations of A and T 2 and each time calculated the error 
after the application of the full swap pulse. In Fig. [S] we 
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FIG. 4: Comparison of error resulting from a full swap op- 
eration and two consecutive square root of swap operations. 
The error is plotted as a function of upper level detuning in 
quantum dot B. The pulse width is optimized for A = 1, 
t\ y 2 = 0.2, [7 = 4, and r = 6. The error for a single square- 
root-of-swap operation is also shown. 

show how the probability of having an error larger than 
the 10~ 4 threshold depends on the energy level accuracy, 
a a- Two cases are considered, namely, plain and limited 
control of the inter-dot coupling constant (err 2 = and 
0.1 T2, respectively). The data indicates that frequent, 
non-correctable errors will happen if an accuracy in A of 
better than 10 percent is not achieved. 



B. Relevance for Real Quantum Dots 

To make a quantitative estimate of the impact of these 
results, let us consider the double QD setup of Jeong and 
coworkers. 25 In their device, each QD holds about 40 elec- 
trons and has a lithographic diameter of 180 nm (we es- 
timate the effective diameter to be around 120 nm, based 
on the device electron density). The charging energy 
and mean level spacing of each QD are approximately 
1.8 meV and 0.4 meV, respectively (thus U/A ps 4.5). If 
we allow for a maximal inter-dot coupling of ix 2 ~ 0.2 A 
(which yields a level broadening of about 0.25 A), we 
find minimal full swap pulse widths of about 100 ps. 
These values match those used in Fig. [SJ For this case, 
switching times of 10 ps would be long enough to oper- 
ate in the adiabatic regime and also provide an efficient 
and well-defined pulse shape. Thus, the combined times 
should allow for 8-10 consecutive full swap gates before 
running into dephasing effects related to orbital degrees 
of freedom^ While these numbers are yet too small for 
large-scale quantum computation, they could be suffi- 
cient for the demonstration of QD spin qubits. Based 
on Fig. [31 we find that accuracies in T2 of about 10% 
would make operations only limited by dephasing, and 
not by fluctuation-induced errors. However, as shown 
in the inset, even for small QDs (typically having small 
U/A ratios), the occurrence of large errors is quite fre- 
quent when level detuning is large. 

An important issue for multi-electron QDs is their 
strong mesoscopic, sample-to-sample fluctuations in en- 
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FIG. 5: Probability of having excessively large full swap errors 
(percentage) as a function of level spacing detuning. The solid 
(dashed) line corresponds to ar 2 /^2 — (0.1) at U = 4 and 
ti = 0.2. The inset shows how the probability varies for a 
fixed width pulse (T w 90), but different charging energies 
U, when there is a large level-position detuning: a a = 0.2, 
or, = (circles) and moderate level-position and coupling 
detunings: <ja = 0.1, ar 2 /^2 = 0.1 (squares). 

ergy level position and wave-function amplitudes. Our 
results so far indicate how big an effect a given change in 
energy or wave-function will produce; now we go further 
and discuss how mesoscopic fluctuations more generally 
affect a collection of qubits. 

In experiments, several electrodes are placed around 
the QD surroundings and their voltages are used to ad- 
just the lateral confining potential, the inter-dot cou- 
pling, and the coupling between QDs and leads. These 
voltages are external parameters that can be used to mit- 
igate the effects of mesoscopic fluctuations by tuning en- 
ergy levels and wave- functions to desired values. Having 
that in mind, our results indicate two different scenarios 
for QD qubit implementations. 

First, if one is willing to characterize each QD pair sep- 
arately and have them operate one by one, mesoscopic 
fluctuations will be irrelevant. It will be possible, with 
a single parameter per QD, say, to isolate and align en- 
ergy levels reasonably well. Errors can be further mini- 
mized by decreasing the inter-dot coupling (thus increas- 
ing T). But since QDs are not microscopically identical, 
each pair of QDs will require a different pulse shape and 
duration. Multi-electron QDs are tunable enough, easy 
to couple, and much easier to fabricate than one-electron 
dots; therefore, multi-electron QDs are most appropriate 
for this case. 

Second, if the goal is to achieve genuine scalability, one 
has to operate qubits in a similar and uniform way, uti- 
lizing a single pulse source. In this case, T and r should 
be the same for all QD pairs. Based on our results above, 
one should strive to maximally separate the top most oc- 
cupied state from all other states, occupied or empty, so 
as to reduce the possibility of leakage during operations 
with a fixed duration. At the same time, it is impor- 
tant to reduce inter-pair cross-talk induced by capacitive 
coupling between electrodes, as well as all inter-dot cou- 
plings except between the top most states of each QD. 
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One should bear in mind that not all electrodes act inde- 
pendently - in most cases a search in a multidimensional 
parameter space has to be carried out. Thus, four tun- 
ing parameters per QD may be necessary to achieve the 
following goals: (i) find isolated, single-occupied energy 
level (two parameters); (ii) align this level with the cor- 
responding level in an adjacent QD (one parameter); (iii) 
control the inter-dot coupling (one parameter). For pa- 
rameters involved in (i) and (ii), an accuracy of a few 
percent will likely be required. Finally, control over the 
inter-dot coupling parameter, (iii), must allow for the ap- 
plication of smooth pulse shapes in the picosecond range. 
Our simulations also show that the pulse width must be 
controlled within at least 0.5% accuracy. Although these 
requirements seem quite stringent, recent experiments in- 
dicate that they could be metiSI 

V. CONCLUSIONS 

In summary, our analysis indicate that mid-size QDs, 
with ten to a few tens of electrons, while not allowing 



for extremely fast gates, are still good candidates for 
spin-qubits. They offer the advantage of being simpler 
to fabricate and manipulate, but at the same time re- 
quire accurate, simultaneous control of several parame- 
ters. Errors related to detuning and sample-to-sample 
fluctuations can be large, but can be kept a secondary 
concern with respect to dephasing effects provided that 
a sufficient number of independent electrodes or tuning 
parameters exists. 
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